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The inspiral and merger of a binary black hole system generally leads to an asymmetric distri- 
bution of emitted radiation, and hence a recoil of the remnant black hole directed opposite to the 
net linear momentum radiated. The recoil velocity is generally largest for comparable mass black 
holes and particular spin configurations, and approaches zero in the extreme mass ratio limit. It is 
generally believed that for extreme mass ratios 77 <C 1, the scaling of the recoil velocity is |V| oc n 2 , 
where the proportionality coefficient depends on the spin of the larger hole and the geometry of 
the system (e.g. orbital inclination). The small recoil velocity is due to cancellations: while the 
fraction of the total binary mass radiated away in gravitational waves is 0(ri), most of this energy 
is emitted during the inspiral phase where the momentum radiated integrates to zero over an orbit. 
Here we show that for low but nonzero inclination prograde orbits and very rapidly spinning large 
holes (spin parameter a* > 0.9678) the inspiralling binary can pass through resonances where the 
orbit-averaged radiation-reaction force is nonzero. These resonance crossings lead to a new contri- 
bution to the kick, |V| oc rf^ 2 . For these configurations and sufficiently extreme mass ratios, this 
resonant recoil is dominant. While it seems doubtful that the resonant recoil will be astrophysically 
significant, its existence suggests caution when extrapolating the results of numerical kick results to 
extreme mass ratios and near-maximal spins. 

PACS numbers: 04.30.Db, 04.25. Nx, 04.70.Bw 



I. INTRODUCTION 

It has long been recognized that gravitational waves 
from an asymmetric source can result in a net emission 
of linear momentum and a consequent recoil or "kick" of 
the system [l|. One of the most important realizations 
of this scenario is in the merger of a binary black hole 
system Q , which has attracted the interest of the astro- 
physics community due to the possibility of disturbing or 
ejecting massive black holes from the centers of galaxies 
and globular clusters 

In recent years, several approaches have been used 
to compute binary black hole kicks. For comparable 
masses and in the strong field part of the merger, numer- 
ical GR computations |13l-l28j now provide the best tool. 
The post-Newtonian series is appropriate for computing 
the waveform emitted during inspiral, and at sufficie ntly 
high order can even follow the plunge phase [1^, 1 3 di | . 
For extreme mass ratio inspirals (EMRIs), one may use 
black hole perturbation theory (BHPT), which treats the 
smaller ( "secondary" ) black hole as a test particle in the 
metric of the primary hole, with slowly varying constants 
of motion to account for radiation reaction followed by 
geodesic motion during the plunge phase [3ll - l34j . One 
may also apply BHPT to the ringdown phase in the form 
of the "close limit approximation," which is based on 
perturbations around the final black hole state [H, [36[ . 
There have also been kick computations [37], HH using the 
"effective one-body" (EOB) method, which maps the fi- 
nite mass ratio problem to the motion of a test particle in 
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a modified black hole metric whose form is constrained by 
agreement with the post-Newtonian expansion |39| : and 
there are computations based on BHPT waveforms but 
with source terms determined from EOB inspiral trajec- 
tories [40( . Most of the analyses have focused on inspiral 
from initially circular orbits, since at large radii gravi- 
tational radiation tends to circularize the orbits on less 
than a merger timescale; but the same techniques have 
been applied to mergers with eccentric initial conditions 

a 

The computational expense of some of these ap- 
proaches has motivated several researchers to propose 
kick "fitting formulae" that return the recoil velocity V 
of the remnant black hole as a function of initial masses, 
spins, and orbital parameters [H, H3, [H, Ii3 - |45j . Ide- 
ally, such a formula would both reproduce numerical GR 
computations and approach BHPT results in the extreme 
mass ratio limit. It should also respect rotation, reflec- 
tion, and interchange symmetries [46l |47| . 

The extreme mass ratio limit of black hole merger kicks 
is perhaps less important from an astrophysical perspec- 
tive because the kicks are small compared to the potential 
wells of galaxies. Nevertheless, it is still of substantial 
theoretical interest: the (moderately) large mass ratio 
case is already starting to enable a comparison of pertur- 
bation theory to full numerical simulations in a regime 
where both are practical and valid (e.g. [H], IHf); and 
BHPT arguments are used to set the limiting behavior 
of the kick fitting formulae, which are often used by as- 
trophysicists. Given that the fitting formulae are typi- 
cally constrained with only a modest suite of simulations 
that sparsely samples parameter space, it is important to 
check their validity in any regime possible. 

BHPT arguments typically give kicks that scale as 
oc rj 2 for EMRIs, where rj = /j,/M <C 1 is the mass ratio 
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[31j . (Since we work to lowest order in /i, we will not 
distinguish here between rj and the traditional mass ra- 
tio q — TOsmaii/wi ar ge; and we will take \x and M to be 
the smaller and larger masses, respectively.) The reason 
is that even though the total energy emitted during the 
inspiral is ~ fj,, emitted in ~ rj^ 1 cycles, in problems 
considered thus far it is emitted symmetrically: averaged 
over an orbit, the amount of power radiated in direction 
n is equal to that radiated in direction — n. resulting 
in zero net recoil. The exception occurs as the inspiral 
ends and terminates in a plunge and finally a ringdown, 
which occurs over only of order one dynamical timescale 
(~ M) and emits an amount of gravitational wave energy 
~ /i 2 /M in an asymmetric pattern. This leads to a kick 
only of order V ~ (/i 2 /M)/M ~ rj 2 , rather than rj, since 
only the energy radiated in the last dynamical time pro- 
duces a unidirectional force. In the case of nonspinning 
black holes as rj — > 0, recent perturbative computations 
have found V = 0.0446r? 2 [H and 0.044V in agree- 
ment with the estimate of 0.04396?y 2 I4Q| obtained via 
extrapolation of numerical simulations [17, 48] . The case 
of EMRIs in the equatorial plane of a spinning primary 
hole also leads to V oc rj 2 , but with a coefficient that 
depends on the primary spin a* [341 ] . 

The principal purpose of this paper is to argue that in 
some cases, an EMRI can produce a kick with a limiting 
behavior V ~ "if I 2 instead of ~ rj 2 . The required condi- 
tions - at least for circular inspirals, which are the sole fo- 
cus of this paper - are (i) a very large spin for the primary 
black hole, a* > 0.9678; and (ii) a low but nonzero in- 
clination prograde orbit. Under such circumstances, the 
inspiralling binary can pass through resonances between 
the vertical and azimuthal frequencies: specifically, the 
increase in longitude between successive ascending node 
passages $ can be an integer multiple of 47r (instead of 
having $ = 2tt as occurs in any spherically symmetric 
spacetime, such as Schwarzschild). When the system is 
in such a resonance, the orbit-averaged recoil force is not 
zero, but rather ~ rj 2 . The EMRI cannot be trapped in 
this resonance since there is no preferred longitude in the 
problem, but rather it continues its inspiral; as a result, 
the resonant argument ip switches its direction of circu- 
lation (ip changes sign). During the resonance crossing, 
the resonant argument has roughly constant phase (vary- 
ing by < 1 radian) for a duration of time t v ~ |</j|~ 1//2 ; 
since tp is a combination of orbital frequencies (of or- 
der M _1 ) that varies on the inspiral timescale M 2 /fi, we 
have (p ~ [i/M 3 and hence t v ~ {i~ 1 / 2 M 3 / 2 . Putting 
this together, we find that the resonant kick is the orbit- 
averaged force times the dephasing time, divided by the 
mass of the system: V ~ (■rj 2 )(fi- 1 / 2 M 3 ^ 2 )/M ~ tj 3 / 2 . In 
this paper, we will use the method of stationary phase 
to demonstrate this scaling, and explicitly evaluate the 
prefactor for some values of spin parameter a* and incli- 
nation L. 

The resonant kick dominates over the 0(j] 2 ) transi- 
tion/plunge kick for sufficiently extreme mass ratios, in 
those cases where it occurs. In fact, for most of the pa- 



rameter space in (a*, l) where a resonance crossing occurs 
it appears likely that the resonant kick will dominate. We 
note that in cases where both kicks are comparable, the 
overall magnitude of the kick will depend very sensitively 
on initial conditions, because - as the vector sum of two 
contributions - it will depend on the relative longitude 
of the ascending node at the resonance crossing and the 
longitude at plunge, as well as the phase of the vertical 
oscillation at plunge. 

None of the currently published kick velocity fitting 
formulae contain an order rf/ 2 contribution. While it is 
not clear whether the resonant kick (or its intermediate 
mass ratio analogue) is significant for astrophysical cases, 
the existence of the if I 2 scaling in some part of (a* , l) 
space suggests caution when constructing fitting func- 
tions or extrapolating numerical GR results to extreme 
mass ratios and/or spins. 

We note that previous work on EMRIs has identified 
the radial-vertical resonances in generic (eccentric and 
inclined) inspirals as potentially important for waveform 
computation: they yield a deviation from adiabatic in- 
spiral behavior at resonance crossings due to radiation 
reaction and self- force corrections [49|, [5(| , or in space- 
times that deviate from the Kerr solution [5l| . However 
these resonances do not exist for circular orbits, and none 
of these analyses appear to have considered the effect on 
the radiated linear momentum. 

This paper is organized as follows. The theoretical 
arguments are presented in SjTTJ The computation and 
the resulting kick magnitudes are given in ^Illi and their 
significance is discussed in i jiVl The formalism and asso- 
ciated code are described at length in Ref. [52J ; we repeat 
only the most important points here and refer the reader 
to Ref. [52j for implementation details. We use relativis- 
tic units where G = c = 1. 



II. THEORY 

We treat the motion of the smaller black hole in the 
Boyer-Lindquist coordinate system 53], in which the 



Hamilton- Jacobi equation for the particle motion [54 1 
and the Teukolsky equation [55l - l58j for the Weyl tensor 
perturbations are separable. The metric in this system 
is 



ds 2 = 



2Mr\ , o h .u in 

dt 2 " : - 2 

E / E 



.2 , „2\2 A„2„:„2 



dt z ^^sm 2 9 dtd<p 



(r 2 + a 2 ) 2 - Aa 2 sin 2 6 , 
sm 2 i 



-— dr 2 + Erf6» 2 , 
A 



2Mr 



and E = r 2 + a 2 cos 2 



(1) 



This 



reduces to standard spherical Minkowski coordinates at 
large r. The mass of the hole is M and its specific angular 
momentum is a < M; we denote the spin as a fraction of 
the maximal value by a+ = a/M. 
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Greek indices a/3... will run over all four coordinates 
{t, r, 9, (f>}, Latin indices ij... will run over only the spa- 
tial coordinates {r, 9, 0}, and the overdot will denote a 
derivative with respect to coordinate time, ' = d/dt. 



A. Geodesic motion 

We treat the motion of a test particle in the Kerr space- 
time using the 3+1 Hamiltonian formalism and utilize a 
canonical transformation from the original {x l ,pi) vari- 
ables to action-angle variables, following the notation and 
methodology of Ref . [|32j . Note that this results in a dif- 
ferent set of angle variables than Refs. [5(J Hll, whose 
4-dimensional analyses use the proper time or Mino time 
[60l . [6l| and promote t to be a dynamical variable with a 
conjugate momentum p tl although the actions Ji are the 
same. 

In this picture, the action is given by the usual test 
particle formula S = J L dt, where the Lagrangian is 
given by L = —^dr/dt, where dr is the proper time ele- 
ment. The action is varied with respect to the trajectory 
{r(t),6(t),4>(t)}. The conjugate momenta {p r ,Pe,P<t>) 
are easily seen to equal the covariant momentum com- 
ponents, and the Hamiltonian H = pii 1 — L is easily 
seen to equal H = —pt, the covariant f-component of 
the momentum determined by the mass-shell condition 
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The Kerr problem admits three constants of the mo- 
tion: the Hamiltonian per unit mass £ — —pt//J-, the an- 
gular momentum per unit mass C = p^/ fJ., and the Carter 
constant Q; and these mutually commute: {£,£}p = 
{£, Q}p = {£, Q}p = 0, where {,}p represents a Pois- 
son bracket. Therefore the particle moves on a 3-torus of 
constant (£, Q,C), which may alternatively be parame- 
terized by the action variables: 



Ji 



1 

2^ 



Ci 



Pj dx° , 



(2) 



where Ci is a loop on a surface of constant (£,Q,C) 
where i € {r, 9, </>} advances through one cycle. One 
may also write the reduced actions Ji = Ji/fJ>, which 
depend only on the trajectory and not on fi. We may 
define the Jacobian of the transformation between the 
two parameterizations of the tori, Mai = 8Ka/ dJi or 
[M _1 ]iyi = dJi/dKA, where Ka is one of £, Q, or C. 
We note that one of the actions is simply the angular 
momentum: = C. 

We also need the angle coordinates on these tori, < 
ijj 1 < 2ir; the mapping from (Jj,ifij) —> (x l ,pi) may be 
explicitly constructed using (i) the direct conditions 



dx l 

dpi 
dip* 



-Aipvr 1 ]^ 



dK A 



dpr 

dK A 



and 



dx l 
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and (ii) an initial condition or choice of origin, i.e. a 
point on the torus to define i\f = ijj e = ip* = 0. There 
is some freedom in choosing the origin, but it is not ar- 
bitrary (the canonical conditions for Ji and tp 1 impose 
some constraints); a valid choice is to set ip' 1 = at peri- 
centre (r = r m i n ), ascending node (9 = tt/2, 9 < 0), and 
zero longitude (0 = 0). 

The Hamiltonian is a function only of the actions and 
so the equations of motion become trivial: J is constant 
and the angles advance at a constant rate, ip l = M$i = 
f2ij we may thus write the angle solution as 



V> = v (0) + sit, 



(4) 



where i/)^ represents the three initial phases. 

In the Keplerian limit, r 3> M, all three frequencies 
become equal, fi r f=a f2e ~ ^<f>- For circular but inclined 
orbits, of interest here, we generally have a precession 
rate (i.e. rate of increase of the longitude of the ascending 
node) Slfjj — fig ^ 0. 

B. Emitted waveform, power, and momentum 

BHPT enables us to compute the waveform emitted by 
an orbiting particle in the Kerr spacetime to linear order 
in its mass fi. It is a practical method of computation for 
EMRIs since it is valid in the strong-field regime, with the 
mass ratio as the only expansion parameter. Computa- 
tion of the orbit-averaged energy and angular momentum 
flux is by this point a "standard" problem and allows one 
to compute the adiabatic evolution of either circular or 
equatorial orbits around Kerr black holes (621-66] . For 
generic (eccentric and inclined) orbits in Kerr, one also 
needs to compute the rate of change of the Carter con- 
stant Q [6(J HI[ , which is one of the more difficult 
problems tackled by BHPT; but since we restrict to cir- 
cular orbits this will not be necessary here. 

The gravitational wave signal is encoded by ip4, the 
perturbation of the Weyl tensor component (not to be 
confused with an angle variable) . It obeys a wave equa- 
tion with a source T(t, r, 6, <p) [561 ] - This equation turns 
out to be separable in the 4 coordinates (the t and <j> 
dependences follow from symmetry arguments, but the 
separation of the r and 9 dependences is nontrivial) , and 
it possesses 3 separation constants {m, uj, A}: those asso- 
ciated with the longitude and time dependences m and u 
(-04 oc e lm * e - JW *) anc i the eigenvalue A of the 9 equation 
(which has boundaries at 9 = and n). For geodesic 
motion, the source is quasiperiodic in the sense that in 
its Fourier transform 7^n,(w,r, 9) contains only frequen- 
cies lj = q • fi, where q € Z 3 is a lattice vector. The 
retarded solution for the gravitational waveform is 



Mt, r, 9, 4>) = (r- ia cos 9)~ 4 



tmq 



xft mw (0)e ! < m ^ <1 > ,O) - ui ». 



(5) 



In this equation: 



4 



S( mu (fl) is an angular eigenfunction normalized by 
Jq \Semui(9)\ 2 sin 8 d6 = 1. It is real, and it reduces 
to a spin- weighted spherical harmonic [iH, [t3| at 
zero longitude, Si m0J (9) =_ 2 Y tm (9,(t> = 0), in the 
limit of a Schwarzschild black hole. The azimuthal 
quantum number m 6 Z, and the vertical quantum 
number £ is an integer with I > max{|m|, 2}. 

7^-3;fanuj( r ) is the purely- outg oing solution to the ra- 
dial Teukolsky equation [5 61 ] ; its large-r behavior is 



(r 



-1 gitjr* 



(6) 



Here r*(r) is a tortoise coordinate whose large-r 
behavior is r* = r — 2M lnr+constant; this ac- 
counts for the logarithmically divergent phase shift 
of emitted radiation in a 1/r potential. 



Resonances 



If / Jig is irrational, then the average value in Eq. (|SJ) 
contains no interference between different (to, fc). If how- 
ever there is a rational ratio Q^/Slg, then different values 
of (m,k) can have the same frequency. We take the res- 
onant relation 



ngflg + n^fi^ = 



(10) 



with ng and relatively prime, so that VL^/VLg — 
—ngjn^. Then there will be interference between two 
terms (mi,fci) and (m^fca) — + Am, k\ 4- Ak) if 
fl^Am + QgAk = 0. This condition is satisfied when and 
only when 



(Am, Ak) = (pn<p 7 png) 



(11) 



• ^m'q i s the emitted wave amplitude associated with 
the Im gravitational wave mode, emitted by the q 
Fourier component (i.e. frequency u) = q • ft). 

We will consider circular ( J r — 0) orbits in our prob- 
lem, in which case there is no dependence of the source T 
on the radial phase; hence we may consider only Fourier 
modes on the torus with q r — 0. Also the longitude 
shift symmetry guarantees that the to Fourier mode of 
the waveform is contributed only by torus Fourier modes 
with (j(fr — TCi. We may thus write q = (q r , qg, q$) = 
(0, k, m). Using this, and the fact that the emitted wave- 
form amplitude is proportional to the particle mass /i, 
we may use the notation = (J-Z^ k , where Z^ k is 

independent of /i. 

[Note that Ref. [52j drops the ip^ term in Eq. §5§ since 
there it amounts to an arbitrary definition of longitude, 
but in our problem we do not have this luxury] 

The power emitted to future null infinity per unit solid 
angle is 



dE 
dtdH 



"^4 dt 



(7) 



where the average value is meant to be taken over many 
cycles. Here J ip4 dt can be obtained by inserting a factor 
of i/u in Eq. (0: 



E 

trnk 



Irak q jo\ i{mij>— q-i/j (0 '— oji) 



dE _ ii 2 
dtd 2 h 4w 

(8) 

The ' in the summation means that we drop terms with 
u = 0, as these do not correspond to any radiated power. 
The emitted power and force are then 



E r . 



dE 9 f dE , 

-— — d n and F = / n , ,„„ d n. (9) 
52 dtd z n J S 2 dtd^n 



The back- reaction force on the inspiralling binary is — F 
in accordance with conservation of momentum. 



for some p £ Z. In this case the two modes have the same 
frequency, 

ui = mifi^, + kiSlg = m 2 fi0 + k 2 flg- (12) 
We may now expand Eq. ([5]) in the resonant case as 



dE 



dtd 2 h 4ir 



— 2 yout* 17' 



E 

£\m\k\£.2P 

x Se imiLJ (0)Se 2rri2UJ (8) 



out 



x exp[ipn - ipn^A^ - ipngip e(0) }. (13) 

In this equation the ^ superscripts on ip l are inconse- 
quential because n^Vt^ + ngflg = and hence the phase 
factor is unchanged if we replace — > 



Air ' 



dtd?n 47r 



— 2 yout* yout 

W Z, £ 1 m 1 fci Z, £ 2 ni2fc2 



£im 1 k 1 £ 2 p 
X Se imi u {d)Si 2m2U j (6) 

x exp(ipn,f,(f) — ipn^ip^ + ipngijj 9 ). (14) 

This is the form of the emitted power that we will use. 
It contains both the torus-averaged contribution to the 
emitted power (p = 0) and the contributions associated 
with the resonant orbit only wrapping around part of 
the torus (p ^ 0). Our principal interest will be in the 
emitted power and force; we see that: 

• The S^-integrated power E cm is not affected by the 
resonant (p ^ 0) terms because the longitude inte- 
gral of exp(ipn0<^) vanishes. 

• Similarly, the z-componcnt of the force F z is not 
affected by the resonance. Moreover, the symmetry 
of the torus across the equator then guarantees that 
there is no time-averaged z component to the force 
in the inspiral phase. 

• The in-plane (x and y) components of the force 
can receive contributions only from terms with 
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pricf, = ±1, i.e. terms in Eq. (fT4|) that have nonzero 
0- integral after multiplying by n x or fij,. This im- 
plies that we need only consider the p = ±1 con- 
tributions to i 7 ^ and F y , and moreover that only 
resonances with \n$\ — 1 can contribute. Moreover, 
under reflection across the equator, F x and F y must 
remain unchanged but ip e is incremented by it (the 
ascending node becomes a descending node) ; there- 
fore there can be no contribution to Eq. (fT4|) from 
terms with png odd. Thus the only resonances that 
contribute to F x and F y have \n^,\ = 1 and ng even. 

The resonance condition fEq. 1 10[) with (ng,n^) is also 
satisfied for (— ng,— n$); so in what follows we simply 
choose the positive sign for n$ = +1. 

The resonances at which there is a net force on the 
binary are thus those where fy : fig are in the ratios 2:1, 
4:1, 6:1, etc., corresponding to ng — —2, —4, —6, etc. 
All of these resonances can in principle occur, but only 
for rapidly spinning primary holes, prograde orbits, and 
small radii. We can see this by considering a maximally 
spinning black hole, a* = 1, and small inclinations i <C 1, 
for which I7li 



fy 
fir, 



1 - 4(r/M)- 3/2 + 3(r/M)~ 



-1/2 



(15) 



this approaches 1 at large radii and oo as r — > nsco = M. 

In practice, evaluating the ratio fy/fy at the ISCO for 
different a* enables us to determine the minimum black 
hole spin at which each resonance is possible. Thus we 
see that the 2:1 resonance can exist for a* > 0.9678; the 
4:1 resonance can exist for a* > 0.99722; and the 6:1 
resonance can exist for a+ > 0.999253. We thus expect 
that the 2:1 resonance will be relevant over the largest 
region of parameter space. We emphasize, however, that 
in order to obtain nonzero emitted amplitudes for k =/= 
the inclination cannot be exactly zero. 

The resonance locations for two values of a* are dis- 
played in Fig. Q] The value a* = 0.998 is a rough upper 
limit to the spin parameter of a black hole spun up by 
accretion due to capture of negative angular momentum 
radiation [72j]; a* = 0.98 shows a less extreme case that 
still possesses a 2:1 resonance. Note that as a* increases, 
additional resonances appear. 



does not remain constant. Rather, as the particle inspi- 
rals and approaches the resonance it circulates at some 
rate (p = ngQg + fy < 0. (Remember that ng is a 
negative even integer.) As it drifts inward, ip increases 
linearly with time and crosses zero at the resonance. 
Thereafter it is positive and the resonant argument be- 
gins to circulate again. The time over which the res- 
onant argument can be treated as roughly constant is 
w |</5| -1 / 2 oc /i -1 / 2 , and there will be an overall recoil 
kick velocity V = — jFdt/M proportional to /i 3 / 2 . 

The stationary phase approximation provides a more 
mathematical statement of the above argument, and is 
the method by which we may obtain the numerical pref- 
actor. Indeed, the reason that circular, equatorial inspi- 
rals into Kerr lead to only a /i 2 contribution to the kick 
is precisely that there is no Fourier component of the 
emitted force [HJ that passes through a stationary phase 
point. 

In order to determine the magnitude of the kick there is 
no need in this problem for a "self-force" calculation: the 
conservative part of the self- force should yield an 0{fi) 
shift in the resonance position, but it is the dissipative 
part (i.e. the part that changes the actions) that deter- 
mines (p. The self-force is required to obtain the direction 
of the kick from initial data at large separations because 
one needs to know the longitude of the ascending node at 
resonance crossing, and hence the phase evolution must 
be accurate to < 1 radian. In the case of inspirals that 
pass through multiple resonances (2:1 and then 4:1, and 
possibly higher) even the magnitude of the total kick de- 
pends on phase evolution because one must know the 
longitude of ascending node advance between successive 
resonance crossings. These issues are beyond the capa- 
bilities of the code in [52j and hence will not be explored 
here. 

We may now apply the stationary phase approximation 
to the p = ±1 contributions to Eq. fT4")) . It is easily seen 
that the two contributions are complex conjugates of each 
other, so 



dE 



dtdH 



, — 2 vout* 



47T " 

P=±l £ 1 m 1 k 1 e 2 

xSi imiU {e)S t „{9)e i{ *-^ 

+c.c, (17) 



D. Kick magnitude 



where mi = m\ + 1 and k2 = k\ + ng. Integrating the x 
and y components of the force using n x — in y = sin 8 e~'^ 
gives a force: 



In order to calculate the magnitude of the kick received 
at the resonance, we need to know how long the force F 
acts. In the limit of an idealized test particle following 
a true geodesic of the Kerr metric, fi — > 0, the mean 
force F acts for an infinite amount of time. However at 
this point, it is necessary to consider radiation reaction 
effects. In particular, the resonant argument 



(f = ngip + ^ 



(16) 



F x — iF v = —e~ iv 
y 2 



i 



— 2 g^out* ^out 
/ j W Zj e 1 m 1 k 1 ^e 2 m 2 k 2 



l\m\k\l 2 



X-St imiU ,(9)St 2m2 oj(6) sin 2 9 d9. 



(18) 



(The complex conjugate term has no contribution here; 
it contributes instead to F x + iF y .) 

Now as we sweep over the resonance, the integral in 
Eq. (|18p varies slowly; the evolution as one moves slightly 
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FIG. 1: The resonance locations in the (r, t)-plane for a* = 0.98 and 0.998. The vertical dotted line on the left in each panel 
shows the horizon, and the dashed curve shows the ISCO location. The solid curves show the resonance locations. Additional 
resonances appear as we increase a*. 



off resonance is dominated by the e itp term. If we 
Taylor-expand tp as 



(19) 



then we have the time integral 



lv dt 




exp I -i(f - i- sgmp ) . (20) 
<p\ \ 4 



(This is the essence of the stationary phase approxima- 
tion.) Using this to integrate Eq. (fl"8|). and using (p = Afj, 
with A independent of fi, we see that the kick velocity is 



|V| = 



2\A\ 



3/2 



E, 2 >70Ut* r^Ollt 



-xSi imiU1 (e)Si 2Tn2UJ (6) sin 2 9d9 



(21) 



Equation (|2ip is the main theoretical result of this pa- 
per: it establishes the existence of a kick during the inspi- 
ral phase of order /^ 3 / 2 under any circumstance that leads 
to a resonant crossing. Nevertheless, we must still com- 
pute the rescaled phase acceleration A at the resonance. 
It is 



A = > — n s — =- . 



(22) 



Since we consider circular orbits, J,, =0; in this case, the 
rate of energy and angular momentum loss are sufficient 
to follow the remaining two "constants" of the motion 



[66l l73j | . Here the rate of change of angular momentum 
is (in the notation of Ref. (H^ ) 



C = J 6 



imk r ' 



(I- 



down 1 2 



'ln,k\ 



(23) 

where a£ mk depends on the separation constants and 
and Zf°^f n are the wave amplitudes emitted to fu- 



£mk 

ture null infinity and into the future horizon (divided 
by fi). We can obtain £ by making the replacement 

m/{2uj^ nk ) — > l/(2cj^ fc ); and then we can obtain Jg using 
the fact that the fundamental frequencies are the partial 
derivatives of the energy with respect to the actions, and 

hence £ = flgJg + ^l^,J^. Since uj m k — k£lg + mSl^, we 

can infer Jg: 



Jo = -M XI 



imk 



2uj j 



r^out \'Z I \ /y down \2 

^imkl a imk\^£mk I 



(24) 



III. RESULTS 

We are now in a position to evaluate the resonant kick 
contribution. We focus on the 2:1 resonance, which oc- 
curs in the largest region of parameter space. The ra- 
diated coefficients Z°^ k and Zf™ k are calculated here 
using the code of Ref. Inclinations reported here 

are defined using the Carter constant as in Ref. [66| : 
tan i = y/Q/C. The inclination of the initial orbit t in ; t is 
not exactly the same as the inclination i res at resonance 
crossing, but the changes in inclination during the inspi- 
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ral are small so we will not integrate the full trajectory 
through the (r, /,)-plane in this paper [73j . 

As a specific example, we consider an inspiral into a 
primary black hole with spin a* = 0.998 and inclination 
i = 20°. The 2:1 resonance crossing is at r = 1.761M, 
with actions Jg = 0.09485M and J = 1.599M, en- 
ergy £ = 0.7659, Carter constant Q = 0.3386M 2 , and 
frequencies fl r = 0.0555M" 1 , fig = 0.1518M- 1 , and 
0^ = 0.3037M -1 . There is no 4:1 or higher resonance 
crossing in this case because one reaches the last stable 
orbit (f2 r = 0) first. 

Considering modes with I < 12, \k\ < 6, we derive a 
force at resonance of |F| = 4.15 x 10~ 4 r) 2 and an emitted 
power of Eem = 6.74 x 10 _2 7y 2 . (The gravitational-wave 
rocket is not very efficient: the asymmetry of the emitted 
power is only \F\/E cm = 0.006.) This leads to an inspiral 
rate of jg = -0.00794^ 2 M and = -0.197?7 2 M, or f = 
— 0.594ry 2 and i — 0.012?7 2 A/ _1 . (The inclination change 
is small compared to the radius change, in accordance 
with previous work [z3i IZ^I-) The second derivative of 
the resonant argument is (p — 0.129?7M _2 , leading to a 
total kick of V = 2.89 x 10~V /2 - This appears to be 
converged with respect to the maximum £ and \k\: the 
kick increases by 2.4% if we go to I < 16, < 8, and by 
a further 0.4% if we go to i < 20, |fc| < 10. 

We have explored the behavior of the resonant kick 
velocity as a function of the primary spin and inclina- 
tion. The results are shown in Figure [2] as a function of 
inclination for a* = 0.98 and 0.998. We can see the ex- 
pected qualitative result that the kick velocity vanishes 
for equatorial orbits (it is oc i? for small C). As the in- 
clination increases, the resonant kick also increases, until 
one approaches the maximum inclination at which the 
resonance lies outside the ISCO (this is 38° for the 2:1 
resonance and a* — 0.998). There the inspiral velocity |f| 
becomes very large and A — > oo. As a consequence, the 
resonant kick, which according to the SPA is proportional 
to |j4| -1 / 2 , drops to zero at the maximum inclination. 

The analysis we have done here assumes that the de- 
phasing time - i.e. the time over which F deviates by 
7r/2 from its the stationary direction, id = ■k 1 / 2 \^>\~ 1 / 2 - 
is much longer than the period of vertical oscillations, Tg] 
using (p = Afx, we see that this condition is equivalent to 

47T 

^ <<? ^MA7f (25) 

(Since A oc M~ 3 and Tg oc M if we change M while keep- 
ing the dimensionless parameters a+ and l fixed, it follows 
that this condition on r] does not depend on M.) For the 
above example of the 2:1 resonance with a* = 0.998 and 
l = 20°, we find r\ a — 0.014. For mass ratios of order 
?7 ~ r) cl -, computing the resonant kick will require a more 
sophisticated analysis, likely one that actually follows the 
inspiral trajectory rather than using the stationary phase 
approximation. 

The kick velocities associated with the 4:1 resonance 
are much smaller: for a* = 0.998, we find a maximum 
resonant kick velocity of ~ 3 x 10 _5 ?7 3 / 2 at inclinations 



Resonant kick velocities 




5 10 15 20 25 30 35 40 
i (degrees) 

FIG. 2: The resonant kick velocity, with the factor of if/ 2 
divided out. We show curves including modes with I < 12 
and \k\ < 6. 



l ~ 15°. While this still becomes larger than the plunge 
kick (oc rj 2 ) in the extreme mass ratio limit rj —> + , the 
4:1 resonant kick velocities are two orders of magnitude 
less than what we obtain from the 2:1 resonance and thus 
are subdominant. 



IV. DISCUSSION 

We have shown that the passage of an extreme mass 
ratio binary black hole through a resonance can lead to 
a fractional-order kick V oc rf/ 2 that dominates over 
the usual nonresonant inspiral + transition + plunge + 
ringdown kick (oc ry 2 ) for small 77. We exhibited this 
effect for the case of inclined circular inspirals into a 
rapidly rotating primary black hole, since the decay of 
eccentricity during the weak-field phase of the inspiral 
makes the circular orbit the most interesting case. How- 
ever, the resonant kick phenomenon should also occur for 
other types of EMRIs, including eccentric inspirals into 
Schwarzschild black holes, eccentric equatorial inspirals 
into Kerr holes, and generic (eccentric and inclined) in- 
spirals into Kerr. The formalism to describe these cases 
would be similar to that used here, but actual computa- 
tions are left to future work. 

The resonant kick described here is only one possible 
resonant interaction in extreme mass ratio binary black 
hole systems. Flanagan & Hinderer [50[ have considered 
resonances between the vertical and radial motions in 
generic orbits, i.e. where n r Q r + ngflg = for n r , ng e Z. 
Unlike the resonances considered here, the phase at which 
one passes through a vertical-radial resonances actually 
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has an influence on the change in the constants of mo- 
tion {£, Q, £), and hence leads to large (oc r\~ x l 2 ) changes 
in the subsequent phase evolution of the inspiral, which 
would be significant for template construction. The res- 
onant recoil effect considered in this paper does not lead 
to such an intrinsic phase shift. It does lead to a Doppler 
shift of the template before versus after the resonance 
passage, but since the fractional frequency shift is oc rf'l 2 
and the number of cycles between resonance passage and 
plunge is oc tj^ 1 , the overall phase shift is oc 77 1 / 2 and 
hence will be small in the EMRI limit. 

It seems unlikely (although not impossible) that the 
resonant kick is directly relevant in any astrophysical 
applications. In the example of Section Mil of primary 
spin a* = 0.998 and inclination 1 = 20°, the kick is 
1 .4(77 /?7cr) 3 ^ 2 kms -1 with rj cr = 0.014 being the maxi- 
mum mass ratio for which we expect our analysis to hold 
(Eq. [23)1 . If rj <C r] cr this is small compared to the veloc- 
ity dispersion of any system that could conceivably host 
a binary black hole. Depending on the actual behavior 
of the resonant kick when rj/r] cr is of order unity, it might 
be significant for globular clusters in the case of a stel- 
lar mass and intermediate mass black hole binary. Even 



then one would have to arrange for a very large primary 
hole spin and the appropriate geometry. 

Nevertheless, the existence of the new rfl 2 kick contri- 
bution suggests that the existing kick fitting formulae be 
treated with caution, particularly in the large-to-extreme 
mass ratio limit. It also highlights the importance of 
exploring the (j < 0.1 regime with numerical GR simu- 
lations [3, [75l477| ; this may reveal the behavior of the 
resonant kick at 77 ~ r/ cr and elucidate the transition from 
rj 2 to 77 s / 2 scaling. 
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